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Abstract 

The numerical algorithm of the inverse quantum scattering is developed. This algorithm is based 
on the Marchenko theory, and includes three steps. The first one is the algebraic Pade approximation 
of the unitary S-matrix, what is realized by solving a system of linear equations. Second step is 
the exact solution of the Marchenko equation. The used approximant reduces it to another system 
of linear equations. At this step we get the real-valued potential. It is shown numerically that the 
developed algorithm is able not only to generate the given S'-matrix dependence, but converges to 
the initial potential. At third step we construct the optical complex-valued potential which gives 
the needed S-matrix. It is shown that the modern phase shift analysis data allow to construct the 
nucleon-nucleon optical potentials of two kinds. These potentials describe the deuteron properties 
and the phase shift analysis data up to 3 GeV and they have different behavior at short distances. 
One is a repulsive core potential and another is a Moscow attractive potential with forbidden states. 



1 Introduction 

Quantum inversion has many applications in nuclear physics. Most of the potential descriptions of the few- 
body nuclear scattering base on the inversion of the scattering data. Though nucleon-nucleon interaction 
is used as input for all nuclear calculations, here we always have some fitted parameters and this fitting is, 
of course inversion. Moreover all modern high-precision nucleon-nucleon potentials are fitted to scattering 
data and are now perceived as phenomenology. The problem here is "that quantitative models for the 
nuclear force have only a poor theoretical background, while theory based models yield only poor results" 
PQ . At the same time the nucleon-nucleon phase shift analysis data are smooth in all investigated energy 
region up to 3 GeV 2 j. This fact justifies the potential description without explicit internal degrees of 
freedom. It is commonly supposed that nucleon-nucleon potential is nonlocal. However local configuration 
space potentials simplify nuclear calculations greatly. Nonlocality effects may be treated as corrections 
due to the internal degrees of freedom. For these reasons construction of the high-precision local NN 
potential describing at least NN scattering data is necessary for exact calculations in nuclear theory. 
Recently such results appeared in the literature [3]. To simplify the following investigations we worked 
out the presented simple algorithm which allows to invert scattering data above inelasticity limit and to 
get the corresponding optical potential. 

In Sect. 2 we describe our inversion algorithm, that allows to get the configuration space potential 
from phase shift analysis data neglecting inelasticity. In Sect. 3 we show how to get an optical complex- 
valued potential from the real-valued one. This optical potential describes phase shift analysis data and 
loss of flux due to inelastic processes. In Sect. 4 we apply the worked out method to the nucleon-nucleon 
scattering data. We extract two different optical potentials corresponding to the different asymptotic 
behavior of the scattering data. 



2 Inversion algorithm 

The input data of the Marchenko inversion are 

{S (q) , (0 < q < +oo) , q j: Mj, j = l,...,n h }, (1) 

where S (q) - is the scattering matrix dependance on the momentum q, q 2 — Em, q 2 = mEj < 0, Ej is 
the energy of the j-th bound state, so that iqj > 0, m - is the particle (reduced) mass. The Mj matrices 
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give the asymptotic behavior of the corresponding normalized bound states. These are output data of 
the partial scattering Schrodinger equation 



d 2 1(1 + 1) \7f X 

ar r 



i>{r, q) = q 2 lb(r, q) (2) 



To illustrate the worked out algorithm first we consider the one channel case. We proceed from the 
Marchenko equation 

+00 

F{x,y)+L(x,y) + J L(x,t)F{t,y)dt = 0, (3) 

X 

where 

, +?° 

F(x,y) = — / h+(qx)(I-S(q))h+( q y)dq + J2 M M(iq j x)h+(iq j y), (4) 

/ij 1 " (z) are the Riccati-Hankel functions. 

Solution of the eq. © is the function L(x,y), which gives the reconstructed potential V (r) for the 
eq. © 

v{r) = _dLpr)__ 
dr 

Contrary to the the algorithm used in j2j we use the algebraic approximant not for the phase shifts 
S(q) but for the S- matrix 

h (q) + ill (q) 



S {q) = ^=( f ^f ) - i ^f ) \ . (7) 



or 

.<? (n\ = P m = ( 

M(q)+ifi(q), 

In the approximant (jSJ fi (q) and ji (q) are an odd and even polynomials of q, which do not turn 
to zero at the real axis simultaneously (for I = 0, in the case then there is a bound state with qi = 0, 
conversely fi (q) and ji (q) are even and odd polynomials). In the approximant Q /1 (q) and ji (q) must 
be of different parity, but we cannot describe the case with zero energy bound state. The important 
feature of approximants © and JJJ) is that these are the most common Pade approximants for the S- 
matrix consistent with its properties. Both approximants lead to the finite-dimensional kernel F (x, y) 
of the equation Q. The function L{x,y) and potential V (r) are expressed through the elementary 
functions (sin(r), cos(r) and powers of r), so we get Bargmann potentials. Choice © or Q) is fixed for 
needed (phcnomenological) dependance S (q) by numerical experiment. We choose the relation which 
approximates S (q) better with less number of S (q) poles (taking into account their multiplicity). The 
less number of these poles gives the more simple potential. 

The properties of the Pade approximants are well known, and there are broad enough functional 
classes, for which this approximation converges everywhere besides poles of these functions. In partic- 
ular, the consequence of the Pade- hypothesis @] is that if a function is analytical at point q = and 
meromorphic at circle D, containing this point, then its diagonal Pade approximations [M/M] (where 
M is the number of the approximation poles) converges to this function on compact subsets of D, which 
do not contain the poles of this function. Obviously for the major part of nuclear physics problems 
the iS-matrix does satisfy these conditions. Particularly for the short range potentials the S'-matrix is 
meromorohic at all complex plane. 

Another remark concerns the S-matrix poles positions. Our approximants for S'-matrix ©, J7J) 
suggest that S'-matrix poles may be positioned everywhere besides the real axis. Previously at numerical 
approaches the following S-matrix approximant was used [S]: 



S(g) = H 



JV 

q 



. q + ia v 
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with real values a v , what assumes that S- matrix poles are positioned only at the imaginary axis. It is 
obvious, our approximants are more common and allow to reconstruct the potentials of more complex 
form. 

The approximants © and Q lead to the following expressions for the phase shifts S (q) 

tg(-Hq)) = Tr\ (») 

or 

»(>>)-£$• 

Having the scattering data set <5 (ft) = ft (i = 1, N), we get system of N linear equations which 
define N unknown coefficients of fi (q)+ifi (q) from J8J) and @. The number of experimental phase shifts 
differs generally from the chosen number N . Therefore the phenomenological dependance S (q) must be 
approximated by appropriate spline. The TV needed values S (ft) = ft are defined by some selection from 
this spline. 

The increase in number of precise experimental phase shifts would make approximants (|SJ|, Q more 
close to the true S-matrix. The scattering data are defined experimentally and they have errors. The 
arbitrary selection of 6 (ft) = ft (i = 1, N) from the experimentally defined allowable region may lead 
to close roots of the polynomials /i (q) and fa (q), i.e. to false S'-matrix poles. Therefore j\ (q) and ji (q) 
must be factored and close roots particularly those situated near the real axis must be eliminated. This 
elimination corresponds to elimination of the Pade approximant defects. 

Furthermore, using the approximant (JSJ (in case of Q the solution is a particular case of the considered 
below solution for two bound channels) and calculating the integral in the eq. 10} using the residue 
theorem, we get the following expression for the kernel F (x, y) of eq. 

F (s, y) = £ Mfhf (ftx) h+ (ftj/) + j2 hht fax) h+ fay) = £ b 3 hf fax) h+ fay), (10) 

i— 1 i—1 J — 1 

where /3j (i — l,...,n po8 ) are the S'-matrix poles. Summing is taken over all the S-matrix poles n pos 
positioned above the real axis and over bound states n h : n = n b + n pos . 

So in this case the kernel of the eq. © is finite-dimensional. It is known that solution of such 
equations reduces to the linear equations solving. We present the specific solution of this equation for 
our case. 

We search the solution of eq. © in the form 

n 

L{x,y)=Y,Pi{x)htfay), (11) 

i=l 

where Pi (x) are unknown coefficients. Substituting i|l(J|) and into © we get the following system of 
linear equations 

n 

Y J A tk {x)P k {x)=D l {x) {i = l,..,n), (12) 
fe=l 

where 



Aik = S ik ~ h ^ ~m ■ > Di{x) = -bihf fax) . (13) 



Afcjt! fax) h+ fax) - (3 k h + fax) h+_ 1 fax) 
The functional coefficients Pi (x) are defined by 112JI 



P t (x)=(A- 1 D) i , (14) 

then L (x, y) and the potential V (r) from (fTTf) and (J5J. 

The worked out technique of the inversion problem was approved on restoration of the square pit 
potentials. The accuracy of the results is illustrated by initial and restored potentials in fig.l. The 
convergence of the restoration algorithm with increasing of power N of the [N/N] Pade approximant of 
the S-matrix is also shown here. 
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In case of two bound channels the system of the partial Schrodinger equations is 

(15) 




where V\ (r), V2 (r) are potentials in channels 1 and 2, Vr (r) is potential bounding them, Xi( r ) an( i X2(?*) 
are channel wave functions. 

By analogy with J7J we approximate the ^-matrix by the following expression 



S(x) 



exp (2iSi) cos2e i exp (i (Si +62)) sin2e 

iexp (i (Si + S2)) sin2e exp (2i^) cos2e 

( f i^m-^m Y (Ji 12) w)'-(^ 12) w)' 2 . jfflwglw n ^hi^ \ 

\ "°(/r > (9)) a +(A (12) (?)) a ,il2 / a 0) (9)+iA W) (9) V^W+i/f'C?)/ (/< 12) (?))V(/ 1 (12) (?)) 2 / 

Here we again choose the most general form of Pade approximant in contrast to the form used in 
The positions of ^-matrix poles are determined from the equations analogous to © , but with additional 
equations for the mixing parameter 

tg(e(q)) = ^\, (18) 

h (q) 

with scattering data e (qi) = E\. 

In case of bound channels and U 7^ 0, in the original Marchenko theory it was proposed to use 
transformation of the initial eqs. Ijl5|l with l\ = 1% = 0. The same approach was used in 0. Nevertheless, 
from the numerical point of view this approach is less effective than the direct solution of the inverse 
problem with generalized Marchenko equation Formally it has the former view 

+00 

L(x,y)+F(x,y) + J L (x, t) F (t, y) dt = 0, (19) 

X 

but functions involved are matrices (2 x 2) 

F(x,y) = — / H(qx)[l-S(q)]H(qy)dq + Y J H(fox)M l H(foy), (20) 

where 

Using the chosen S'-matrix approximant and the residue theorem we get the following expression for 
F(x,y) 

F(x,y)=i £ ResH (qx) (I - S (q)) H (qy) + J2 H (fox) (foy) = 

Imf3i>0 i=l 

= £ H(fox)Q\H(foy) + £ xH' (fox) CftH (foy)+ (22) 

Jm/3;>0 ftGA 

+ £ H(fox)Q^H'(foy)y, 
fteA 

here (j = 1, 2) are constant matrices, 
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[ > ~ \ dh+(x)/dx 



Pi are the .S'-matrix poles, values qj correspond to bound states, A is the set of S'-matrix poles of the 
second order positioned above the real axis. 
We solve eq. (|19|l using substitution 

n n 

L (x, y) = ^P ( X ) H (PiV) + J2 N * ^y H ' ' ( 23 ) 

»=1 i=l 

where summing is taken over all the S'-matrix poles n poa positioned above the real axis and over all values 
qj. This construction leads to the system of linear equations for the functional (2 x 2) matrix-coefficients 

P i (x),N i (x) 



E Pi ( x ) Qij ( x ) +T,Ni (x)Q% (x) = H fax) Q) + xH' (p jX ) Q) 

E k 0=) Qlj 0=) + E Pi {x)Q% (») = h (p jX ) Q) 



(24) 



where 



(x) = ISij + J H (Pit) H (fat) dt x Q) + J tH (fat) H' (fat) dt x Q 2 

X X 

Q%(x) = [ H (fat) H (pjt) dt x Q 2 (25) 



+00 



tH' (Pit) H (fijt) dtxQ)+ / t 2 H' (fat) H' (Pjt) dt x Q 2 



Q% (x) = ISij + J tH' (p l t)H(P 3 t)dtxQ 2 , 

X 

I is the unit matrix. The integrals in expressions (|25J) are easily calculated analytically, but are cumber- 
some so we do not present them. 

Having solved this linear equation system we get the sought-for potential from (12;it and JSJ). 

The multichannel generalization is made analogously. 



3 The optical potential 

The Marchenko inversion does not give the needed optical potential. But found potential may serve as 
the initial potential for the iteration procedure that converts it into the optical (complex- valued) potential 
that describes inelastic processes. 

First we consider the one channel problem. 

The phase equation [7] for the initial potential V° (r) obtained by some inversion procedure (say from 
Marchenko equation) is 

oo 

5f } = -\ f (r) Df (qr) sin 2 (qr + 6 (0) (r)) dr, (26) 
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where Di (z) is Riccati-Bessel amplitude 0. 

Let us consider the complex-valued potential V^ 1 -* (r) obtained from V° (r) by transformation 

V w (r) = (1 + ia) T/ (0) (r) , (27) 
where a is some real parameter. Evidently the phase equation for this potential is 



oo 

<5 (1) = ~ (1 + ia) y y (0) (r) D? (gr) sin 2 (gr + 5 {1) (r)) 



dr. (28) 



From eqs. ((2611 and l|28|l we get 

- (1 + ia) = 

= _ I±M J y(o) (r) D 2 (At) (sin 2 (gr + (r)) - sin 2 (gr + (r))) dr = 
o 

= _i^a jV<°> (r) Df (At) sin (2gr + <5« (r) + 5^ (r)) sin (<5« (r) - *(°) (r)) dr 



(29) 



For smooth enough potentials and a = a (q) not rapidly increasing with increasing of q the right side 
of eq. (|29|) rapidly decreases comparing with S^ ' and 5^, because under the integral in 129fl there is 
more rapidly oscillating function than in (|26|l and l)28|) . Then as a first approximation we may take 

<5 (1) » (l + ia)*<°> =5 R + i5 I . (30) 

For inelastic scattering the S'-matrix is expressed through the real inelastic parameter p and the real 
phase shift S 

S = cos 2 (p) e 2lS = e 2i(5lt+iSl) , (31) 

so it is easily arrived at 

6^5(°\ (32) 



cos 2 (p) ^ e- 2aS( °\ (33) 

whence it follows that a5 > 0. 

The formula (|33J) allows to calculate the parameter a from the known values p and ps 5. 

This approximation works well enough for many quantum mechanical problems. For example in case 
of nucleon-nucleon scattering p = below the inelasticity limit which is high enough and p grows slowly 
enough above this limit. 

Nevertheless in some instances the initial potential V 1 - ^ (r) must be corrected. In this case the 
following iteration procedure may be efficient. First from the formula i|33|) values a(q) are calculated 
from p(q) and 5(q). Then with the potential (r), defined by formula (|27|) new phase shifts 5^ are 
calculated from J2HJ>- If they do not satisfy phase shift analysis data 5^ ± AS then changing the input 
data to S® — 2S^ — <jW we repeat inversion (eqs. l(T|l-l|14|l) and determine new a from formula 1)33( 1 . 

The case of two channels is considered in a like manner, though the final expressions are more complex. 
By analogy with the one channel case the following generalization for the optical potential is derived (the 
mixing parameter e is considered small): 

v <i) (P) = ( (1 - *«) < (! - ^ + 0/2) \ 

\ (l-i(a + /3)/2)v£> {l-mvg> j 

The first S'-matrix approximation is 
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g(X) _ f e^ 2l&1 ^> cos2ecos 2 p x e i(Si+5 2 ) s i n 2e cos cosp 2 

\ e l ( <5 i+' 52 ) sin2e( 1 ^ cos pi cosp2 e^ 2l& ^ cos 2e cos 2 pi 

e 2»(l +M )C CQs2 ^ + £ (0) e ,( 5 < 1 '(l + ») + 4 1) (l + ^)) sin2 ^ + .(o+^ £ (0) 

e i^(l+ia)+^(l + i0)) sin2 ^ + .Jo4^ £ (0) e M(l+i0#> CQs2 ^ + iic+aj £ (0) 

(35) 

As all the parameters and phase shifts are real, then eq. I|35|) links between parameters a and (5 with 
the mixing parameters pi and phase shifts Si = 5^ . 

This consideration may be generalized to the multichannel case though calculations become much 
more complex. 



4 The optical nucleon-nucleon potential 



We applied the described algorithm of inversion to reconstruction of the nucleon-nucleon potential. 
As input data for this reconstruction we used modern phase shift analysis data up to 1100 MeV for 
3 Si — 3 Di state and up to 3 GeV for 1 S , state of nucleon-nucleon system The deuteron properties 
were taken from [Sj. 

The relativistic effects were taken into account in frames of relativistic quantum mechanics of systems 
with a fixed number of particles (point form dynamics). The review of this approach can be found in [§]• 

This approach is based on the assumption that at not high energies we may consider the number of 
particles fixed, but the invariance group is the Poincare group. A system of two particles is described by 
the wave function, which is an eigenfunction of the mass operator or the mass squared operator. In this 
case we may represent this wave function as a product of the external and internal wave functions [TXillTT] . 
The internal wave function is also an eigenfunction of the mass operator or the mass squared operator. 
It is shown that the mass squared method is consistent with conventional fitting of the Lorentz invariant 
cross section as a function of laboratory energy [13] . We consider system of two particles (nucleons) with 
equal masses. Then the internal wave function xcQ) satisfies the following equation 

[4(m 2 + q 2 ) + V] X = M 2 X (36) 

or 

(- + T-)x = E X , (37) 
\ m 4771 / 

where 

_ M 2 - 4m 2 _ k 2 
4m m 

The eq. (I3TI) formally coincides with the Schrodinger equation. In eqs. lpT7l) M 2 is the mass 

squared operator, m is the mass of nucleon, V is the nucleon-nucleon potential, q is the momentum 
operator of one of the nucleons in the center of masses system. We use system % = c = 1. The 
quasicoordinate representation corresponds to the realization q = — i^, V — V(r). In 12 we showed 
that this formalism can be easily generalized to the case of inelastic channels, particularly it allows to 
take into account the isobar channels in NN scattering. 

This formal coincidence allows us to apply our inversion algorithm. The modern phase shift analysis 
data P] allow to construct nucleon-nucleon potentials of two different kinds, depending on the asymptotic 
behavior of the phase shifts above investigated energy region of 3 GeV. We constructed the nucleon- 
nucleon optical potentials of two kinds for 1 Sq and 3 Si — 3 D\ partial waves. These potentials describe 
the deuteron properties and the phase shift analysis data up to 3 GeV and they have different behavior 
at short distances. One is a repulsive core potential and another is a Moscow attractive potential with 
forbidden states. The Moscow potential was introduced in These potentials are not phase equivalent 
and their S'-matrices differ even bellow 3 GeV, but within the experimental errors. The 1 Sq phase shift 
of Moscow potential begins from it. 3 Sq phase shifts of Moscow potential begin from 27r. The mixing 
parameter ei of Moscow potential differs from this of repulsive core potential by sign. Above 3 GeV these 
potentials give different S'-matrices. 
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The real parts of the constructed partial potentials are presented in fig. 2 and fig. 3. They may 
be downloaded from www.physics.khstu.ru in numerical form. Standard notation for central and tensor 
parts is used, so for 3 Si — 3 D\ 

Vcs(r) = Vt(r), V CD {r) = V 2 {r) + V T (r)/ 2y/2, V tens (r) = V T {r)/V2 

The imaginary parts of potentials are defined by eqs. 1)271 134|) . where parameters a and j3 can be easily 
calculated from the phase shift analysis data .2 . The phase shifts and mixing parameter are compared 
with the phase shift analysis data [5] in fig. 4. In table 1 the results of our deuteron properties calculations 
are compared with the experimental data [5] . Both kinds of the constructed partial potentials describe 
the data well in the limits of the experimental errors. 

We conclude that the available data of the NN scattering and the deuteron properties are not dis- 
criminative with respect to the kind of the NN potential. To discriminate between the considered two 
kinds of nucleon-nucleon potential we need careful experimental examination of other inelastic processes 
such as pp — > ppj or 2 H + 7 — > n + p. 
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Table 1. The deuteron properties 





Exp. a 


Calculation 
with repulsive 
core potential 


Calculation 
with Moscow 
potential 


Energy (MeV) 


2,22458900(22) 


2,2246° 


2,2246° 


Q (Fm 2 ) 


0,2859(3) 


0,2639 c 


0,277 c 


A s (Fm" 1 / 2 ) 


0,8802(20) 


0,8802 


0,8802 


r d (Fm) 


1,9627(38) 


1,951 


1,956 


Vd/s 


0,02714 


0,02714 


0,02714 




0,857406(1) 


0,8497 c 


0,859 c 



a relativistic correction included; b Data are from 8 ; c Meson exchange currents are not included. 
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Figure 1: Reconstructed square pit potential. Solid line N=31, dotted line N=23. N is power of the 
[N/N] Pade approximation. 




Figure 2: Real parts of optic repulsive core potentials. Solid line 1 So -one channel. Two bound channels 
long dashed line Vcs(r) ( 3 Si), short dashed line Vcd{t) ( 3 Di), dotted line Vt ens (r). 
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Figure 3: Same as in in Fig. 2 but for the Moscow potentials. 
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Figure 4: Phase shifts and mixing parameter. Solid line for reconstructed Moscow potential, dashed lines 
for reconstructed repulsive core potential. The phase shift analysis data are from [5]. For 5* waves there 
are two sets of the phase shifts. First one is the original data set from - small symbols. These data 
are described by the repulsive core potential. Second one are the same phase shifts raised 180 degrees up 
- big symbols. These data are described by the Moscow potential. To leave the ^-matrix unchanged we 
must then change the sign of the mixing parameter ei for the Moscow potential. The mixing parameters 
for both kinds of potentials differ only by sign in our calculation and corresponding curves coincide in 
this figure. 
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